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ABSTRACT 

We have performed hydrodynamical simulations to investigate the effects of galactic winds on the high-redshift 
(z = 3) universe. Strong winds suppress the formation of low-mass galaxies significantly, and the metals carried by 
them produce C IV absorption lines with properties in reasonable agreement with observations. The winds have 
little effect on the statistics of the H I-absorption lines, because the hot gas bubbles blown by the winds fill only 
a small fraction of the volume and because they tend to escape into the voids, thereby leaving the filaments that 
produce these lines intact. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: formation — intergalactic 
medium — quasars: absorption lines 



1. INTRODUCTION 

Feedback from star formation is thought to play an impor- 
tant role in the formation of galaxies. For example, theoretical 
models require feedback in order not to overproduce the num- 
ber of low mass galaxies (White & Rees 1978) and the fraction 
of baryons that cools (e.g., Balogh et al. 2001). Models with- 
out feedback also predict disk galaxies that are too small (e.g., 
Navarro & Steinmetz 1997) and an X-ray background that is 
too strong (e.g., Pen 1999). 

Observations of galactic winds indicate that strong feedback 
processes do indeed occur. Observations in X-rays, and of opti- 
cal and UV lines, indicate that most local starbursts (e.g., Heck- 
man 2000), as well as high-redshift Lyman Break galaxies (e.g., 
Pettini et al. 2001), drive winds with a mass loss rate compa- 
rable to their star formation rate. The wind speeds are high, 
100-1000 km s" 1 , and so the conversion of supernova (SN) en- 
ergy into kinetic energy must be quite efficient. This process 
is currently not well understood because of the complications 
that arise from the multiphase nature of the interstellar medium 
(ISM, e.g., McKee & Ostriker 1977; Efstathiou 2000). Simula- 
tions that take some of these complications into account, have 
shown that SNe can indeed plausibly power a wind (e.g., Mac 
Low & Ferrara 1999). 

Although successful models of galaxies appear to require 
feedback, the same is not true for models of the high-redshift 
(z p> 2) intergalactic medium (IGM). The IGM can be studied 
in great detail from the properties of the hydrogen absorption 
lines, seen in the spectra of background quasars (see Rauch 
1998 for a review). Hydrodynamical simulations of the IGM 
(see Efstathiou, Schaye & Theuns 2000 for a recent review) as 
well as semi-analytic (Bi & Davidsen 1997; Viel et al. 2002) 
and analytic (Schaye 2001) models, have been very success- 
ful in reproducing the statistical properties of the observed H I 
lines. These models, which generally do not take feedback into 
account (but see Cen & Ostriker 1999), suggest that the Lya 



forest absorption arises in a network of voids and filaments, 
with the higher column density absorbers located at the inter- 
sections of filaments. The low column density absorbers are ex- 
tended structures with densities around the cosmic mean, which 
contain a large fraction of the baryons in the universe. 

Although feedback has so far largely been ignored in mod- 
els of the IGM, the detection of metals in the IGM suggests 
that it does play a role. The higher column density (TVhi ;> 
10 14 5 cm -2 ) Lya absorption systems generally have detectable 
absorption by C IV (Cowie et al. 1995) and, at least at z ^ 2.5, 
O VI (Carswell, Schaye, & Kim 2002). Furthermore, there is 
statistical evidence that metals are also present at somewhat 
lower densities (Cowie & Songaila 1998; Ellison et al. 2000; 
Schaye et al. 2000a). Simple photo-ionization models indicate 
that the absorbers have a metallicity of order 0.1-1 per cent 
solar (e.g., Cowie et al. 1995; Rauch, Haehnelt, & Steinmetz 
1997; Hellsten et al. 1997; Carswell et al. 2002). 

If galactic winds are ubiquitous and able to transport mass 
to large distances, then they may be responsible for enriching 
the IGM with metals. Indeed, numerical simulations of galac- 
tic winds in a cosmological setting (e.g., Gnedin 1998; Cen & 
Ostriker 1999; Aguirre et al. 2001a, 2001b; Thacker, Scanna- 
pieco, & Davis 2002; Springel & Hernquist 2002) suggest that 
winds could enrich a substantial volume fraction of the IGM to 
the inferred levels. 

Therefore, both observational and theoretical considerations 
suggest that some fraction of galaxies may undergo an episode 
in which they blow a strong wind into their surroundings. This 
may have observational effects on the Lya forest in quasar 
spectra (Theuns, Mo, & Schaye 2000; Croft et al. 2002), which 
may have already been detected (Rauch et al. 2001; Adelberger 
et al. 2002). In this Letter, we use hydrodynamical simulations 
to investigate whether galactic winds that strongly influence the 
properties of small galaxies, and are effective in enriching the 
IGM with metals, can do so without undermining the success 
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FIG. 1 . — Top panels: Density (a, b) and temperature (c, d) for a slice through the simulation with feedback (a, c) and the simulation without feedback (b, d), at 
redshift 2 = 3. The simulation box is 5.u7i~' (co-moving) Mpc. The slice is chosen to go through the most massive object in the simulation. In the feedback simu- 
lation, hot bubbles of gas surround the galaxies. These tend to expand into the voids, thereby leaving the filamentary network of higher density regions unaffected. 
The four numbered horizontal lines (from top to bottom) are sight lines for which the hya spectrum, temperature, and density, are shown in the bottom panel (from 
left to right, offset for clarity). The feedback and no feedback cases are shown as blue and red lines, respectively. Only the strongest lines are significantly affected 
by feedback. 



of current models of the IGM. 

2. HYDRODYNAMICAL SIMULATIONS 

We have performed simulations of vacuum energy domi- 
nated, cosmologically flat cold dark matter models, with cos- 
mological parameters (fi m ,^/,/! 2 ,/!,fJ8,F)=(0.3,0.019,0.65,0.7,0.24). 
The lower value of erg was chosen in view of the recent re- 
sults from the Two-Degree Field survey (Lahav et al. 2002). 
The simulations are performed using a modified version of 
HYDRA (Couchman, Thomas & Pearce 1995; Theuns et al. 
1998), in which the gas is photo-heated by an imposed uniform 
UV-background which evolves with redshift and reionizes H I 
and He I at z ~ 6 and He II at z ~ 3.2. The photoheating rates 
were adjusted in order to match the temperature measurements 
of Schaye et al. (2000b), where the current simulation is re- 
ferred to as the 'designer' simulation. Non-equilibrium rates 
for photoionization, cooling, and heating are computed using 
the fits in Theuns et al. (1998). The particle masses for dark 
matter and baryons are 6.5 and 1.1 x 10 6 M respectively, and 
the simulation box (which has periodic boundary conditions) is 
5.0h~ l (co-moving) Mpc on a side. The Plummer softening is 
5/T 1 (co-moving) kpc. 

Star formation and feedback are implemented using the pre- 
scription described in detail in Kay et al. (2002). Briefly, cold 
gas with density p/(p) > 1000 and temperature T < 50,000 K 
is converted into collisionless stars. Each star particle heats the 
nearest gas particle using equation (3) of Kay et al. (2002), 
with a given efficiency parameter (we use e = 1). Reheated gas 
is prohibited from cooling for 10 7 yrs to crudely account for 
the multi-phase nature of the ISM. This simple feedback pre- 
scription has the desired effect of blowing hot bubbles around 



star-forming regions. The heated gas is also enriched with met- 
als. A yield of three times solar matches the metal line data 
well (see below). 
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FIG. 2. — Stellar mass function of galaxies at redshift z = 3, for the simu- 
lation with feedback (full line) and without feedback (dashed line). The feed- 
back scheme has a dramatic effect on the masses of the galaxies. 



We make no attempt to model the transition from SN bubble to 
galactic wind correctly, since our simulations lack both in res- 
olution and input physics. We do not include cooling due to 
metal transitions which could be important. Our primary aim 
in this preliminary investigation, is to test whether strong feed- 
back in the form of galactic winds can coexist with the Lya 
forest as we know it from simulations without feedback. We 
use the galaxy mass function and the properties of the CIV 
forest to demonstrate that the implemented feedback is indeed 
significant. We have performed two simulations with identical 
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FIG. 3. — Column density distribution functions (CDDFs) of H I (panel a), C IV (panel b), and the C IV versus H I column density of systems (panel c). Filled 
circles refer to the combined line-lists of quasars Q1422+231 (Ellison et al. 2000; z em = 3.6) and APM 08279+5255 (z em = 3.91), and the solid line indicates the 
results for the feedback simulation at z = 3 . Panel (a) : The H I CDDF of the feedback simulation is nearly identical to that of the simulation without feedback (dotted 
line and triangles), and both fit the observations very well for logNni < 15. Panel (b): The feedback simulation matches the observed C IV CDDF, but imposing a 

uniform metallicity (Z = 1O~ 3 Z0, hashed region) does not work so well. Panel (c): C IV versus H I column density, for systems identified on a smoothing scale of 
± 150km s -1 . The median for the feedback simulation is indicated with the full line, with dashed lines indicating the 5 and 95 per centiles. The hashed region is the 
corresponding range for the uniform metallicity case. The feedback simulation reproduces the observed median and scatter better than the uniform Z case. 



initial conditions, one with and one without feedback, so that 
we can examine its effect directly. 

3. RESULTS 
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FIG. 4. — Line- widths distribution from the combined line-lists of quasars 
Q1422+231 (Ellison et al. 2000; Zem = 3.6) and APM 08279+5255 (z em = 3.91) 
for C IV (filled circles) and H I (open circles). Full lines are the corresponding 
results for the feedback simulation. Dashed lines are the H I and C IV line- 
width distributions for the simulation without feedback (and uniform enrich- 
ment). (Only lines for which 12 < log^Afciv < 13.2, and errors in Nciv and 
bciv less than 0.3 dex are included.) The feedback simulation matches both 
distributions well, and feedback has little influence on the H I distribution. 

Figure 1 illustrates the effect of feedback from star formation 
on the density and temperature of the IGM. In the simulation 
without feedback, the density structure displays the usual fila- 
mentary pattern of mildly overdense regions, mostly at temper- 
atures ~ 10 4 K, although accretion shocks have also generated 
some significantly hotter gas. In the simulation with feedback, 
galaxies are surrounded by hot, windblown bubbles; in particu- 
lar, the most massive object in the slice is embedded in a halo 

7 The discrepancy between the simulations and the observations for Vm > 10 15 5 
the simulation box. 



of ~ 1 Mpc co-moving size, of temperature up to 10 7 K. 

Except for the very strong absorption lines, the Lya spectra 
are not significantly affected by these hot bubbles (Fig. 1, bot- 
tom panel). In contrast, feedback has a very strong effect on the 
stellar masses of the galaxies (Fig. 2). The hot bubbles disrupt 
or prevent inflow of cold gas, which reduces the star-formation 
rate dramatically. 

We have created absorption spectra along random sight lines 
through the simulations at z = 3 and fitted them with Voigt 
profiles using VPFIT (Webb 1987), mimicking a high signal- 
to-noise, high resolution spectrum (see Theuns, Schaye & 
Haehnelt 2000 for details of the procedure). We imposed a 
power-law ionizing background J u oc v a with spectral slope 
a = -1.5 and a break at the He II Lyman limit (4 Ryd) such 
that the softness parameter (i.e., the ratio of the H I and He II 
ionization rates) THi/THen = 10 3 . Such a soft spectrum may be 
appropriate for z <; 3, when He II reionization was incomplete 
(e.g., Theuns et al. 2002 and references therein). The ampli- 
tude of the ionizing flux is scaled to give the simulations the 
observed mean Lya absorption. The required flux differs by 
around 10 per cent between the two simulations, with the sim- 
ulation with feedback requiring more ionizing photons because 
the density in voids is enhanced by the winds. (Note that this 
ionizing background differs from the one imposed during the 
simulation, see Theuns et al. 1998 for tests of such rescaling.) 

The Lya-column density distribution functions (CDDFs; the 
number of absorption lines per unit column density, per unit 
redshift) for the two simulations are compared in Fig. 3a: they 
are nearly indistinguishable and fit the observations well 7 . The 
same is true for the distribution of absorption line widths (b- 
parameters; Fig. 4). This is the main result of this Letter, even 
though the galaxies in the feedback simulation drive strong 
winds, there is no discernible effect on the Lya-forest. Closer 
examination of the properties of the hot bubbles reveals why 
this is so: the winds expand preferentially into the lower den- 
cm -2 is caused mainly by the missing large-scale power due to the finite size of 
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sity regions (Fig. 1), and so keep the filaments that produce the 
hydrogen lines intact. Only very strong lines differ noticeably 
between the two simulations, but this does not affect the statis- 
tics significantly. 

Another way to demonstrate that our feedback implementa- 
tion does influence the IGM is by examining the C IV forest. 
The hot gas surrounding the galaxies has been enriched with 
metals and produces C IV absorption. We used CLOUDY 
(version 94; Ferland 2000) to solve the ionization balance 
equations, computed mock C IV spectra (see Aguirre, Schaye 
& Theuns 2002 for details) and analyzed them with VPFIT. 
Fig. 3b shows that the feedback simulation reproduces the ob- 
served C IV CDDF very well. In contrast, the low and high end 
of the CDDF cannot be matched simultaneously by imposing a 
uniform metallicity (hashed region). 

The metal distribution is very inhomogeneous in the simu- 
lation, since metals are deposited by the winds predominantly 
in the surroundings of galaxies. This causes a scatter in the 
C IV column density for a given Lya column density. We have 
proceeded as follows to quantify this effect. Find the strongest 
Lya-line in the line-list, and sum up the column-densities of all 
H I lines, as determined by VPFIT, of which the center falls 
within ±Akms _1 . This is the HI column density of that 'sys- 
tem' . Use the same algorithm to compute the C IV column- 
density at the corresponding wavelength. Remove all lines that 
have contributed to the column-density from the line list, and 
start again. Fig. 3c plots CIV versus HI column density for 
A = 150. The feedback simulation reproduces the median and 
scatter in the observations relatively well, whereas the uniform 
metallicity case underpredicts the scatter (see also Dave et al 
1998). 

However, some caution is appropriate. Given that in the ob- 
servations Ncw/Nm ~ 10~ 2 6 , with no evidence for a trend with 
Nhi (e.g., Ellison et al. 2000), the lack of strong C IV lines for 
the uniform metallicity case may be caused in part by the fact 
that the simulation underpredicts the number of strong H I lines. 
Using a harder UV background, would increase the metallicity 
required to match the observations. On the other hand, includ- 
ing metal cooling would increase the C IV fraction (and thus 
strengthen the C IV lines) considerably for the feedback simu- 
lation. With our choice of parameters, the feedback simulation 



reproduces the C IV data well. However, other choices of pa- 
rameters might do equally well. 

In the feedback simulation, most of the intergalactic met- 
als reside in hot (T ~ 10 4 5 - 10 6 5 K), relatively metal-rich 
(Z ~ 0.1 — 1.0 Zq) gas, and collisional ionization has a signif- 
icant effect on the strength of metal lines such as CIV, NV, 
and O VI. Consequently, the temperature of the gas that domi- 
nates the absorption in one transition, e.g. C IV, can be different 
from that of another transition, e.g. H I or O VI. In such a sce- 
nario one would expect the metal line widths (or ^-parameters) 
to be higher than in the simulation without feedback (and uni- 
form enrichment). This is in fact not obvious from Fig. 4, be- 
cause the lines also have a significant nonthermal component. 
At present, both distributions seem to be in reasonable agree- 
ment with the data (filled circles with error bars). Again, some 
caution is appropriate: metal cooling (which was not included) 
could reduce both the line widths and the importance of colli- 
sional ionization in the feedback simulation. A more detailed 
investigation of the metal lines is beyond the scope of this letter. 
Here we merely conclude that the feedback simulation predicts 
a (highly inhomogeneous) metal distribution that appears to be 
in reasonable agreement with the current data, but differs sig- 
nificantly from the uniform metallicity models considered pre- 
viously (e.g. Rauch et al. 1997). 

In summary, galactic winds that are strong enough to 
strongly suppress galaxy formation and to pollute the IGM with 
enough metals to reproduce the data, do not have a significant 
effect on the Lya-forest. The reason is that the winds fill only a 
small fraction of the volume and tend to expand into the voids, 
leaving the filaments that produce the hydrogen lines intact. 

Acknowledgments TT thanks PPARC for the award of an Ad- 
vanced Fellowship, and M Pettini and M Rauch for comments 
on the ms. JS acknowledges support from the W.M. Keck Foun- 
dation and MV from an EARA Marie Curie Fellowship un- 
der contract HPMT-CT 2000-00132. This work was supported 
by the European Community Research and Training Network 
'The Physics of the Intergalactic Medium' and was conducted 
in cooperation with Silicon Graphics/Cray Research utilizing 
the Origin 2000 super computer at the Department of Applied 
Mathematics and Theoretical Physics in Cambridge. 



REFERENCES 



Adelberger, K., et al. 2002, preprint 

Aguirre, A., Hernquist, L., Schaye, J., Weinberg, D. H., Katz, N., & Gardner, J. 

2001a, ApJ, 560, 599 
Aguirre, A., Hernquist. L., Schaye, J., Katz, N., Weinberg, D. H., & Gardner, J. 

2001b, ApJ, 561,521 
Aguirre, A., Schaye, J., & Theuns, T. 2002, ApJ, 576, 1 

Balogh, M. L., Pearce, F. R.. Bower, R. G., & Kay, S. T. 2001, MNRAS, 326, 
1228 

Bi, H. & Davidsen, A. F. 1997, ApJ, 479, 523 

Carswell, R. F, Schaye, J., Kim, T.-S. 2002, ApJ, in press ( astro-ph/020437CI) 

Cen, R. & Ostriker, J. P. 1999, ApJ, 519, L109 

Couchman, H. M. P., Thomas, P. A., & Pearce, F. R. 1995, ApJ, 452, 797 
Cowie, L. L., Songaila, A., Kim, T, & Hu, E. M. 1995, AJ, 109, 1522 
Cowie, L. L. & Songaila, A. 1998, Nature, 394, 44 

Croft, R. A. C„ Hernquist, L, Springel, V., Westover, M., White, M. 2002, ApJ, 



submitted ( astro-ph/020446C 
Dave, R., Hellsten, U., Hernquist, L., Katz, N., & Weinberg, D. H. 1998, ApJ, 
509,661 

Efstathiou, G. 2000, MNRAS, 317, 697 

Efstathiou, G., Schaye, J., & Theuns, T. 2000, Royal Society of London 

Philosophical Transactions Series, 358, 2049 
Ellison, S. L., Songaila, A., Schaye, J., & Pettini, M. 2000, AJ, 120, 1175 
Ferland, G. J. 2000, Rev. Mexicana Astron. Astrofis. Ser. Conf., 9, 153 
Gnedin, N. Y. 1998. MNRAS, 299, 392 

Heckman T., 2000, in ASP Conf. Ser. 240, Gas and Galaxy Evolution, series, 
ed. JE Hibbard, MP Rupen & JH van Gorkom (San Fransico: ASP), 345 



Hellsten, U., Dave, R., Hernquist, L., Weinberg, D. H., & Katz, N. 1997, ApJ, 
487,482 

Kay, S. T., Pearce, F. R., Frenk, C. S., & Jenkins, A. 2002, MNRAS, 330, 113 

Lahav O., et al., 2002, MNRAS, 333, 961 

Mac Low, M. & Ferrara, A. 1999, ApJ, 513, 142 

McKee, C. F. & Ostriker, J. P. 1977, ApJ, 218, 148 

Navarro, J. F. & Steinmetz, M. 1997, ApJ, 478, 13 

Pen, U. 1999, ApJ, 510, LI 

Pettini, M., Shapley, A. E., Steidel, C. C, Cuby, J., Dickinson, M., Moorwood, 

A. F. M., Adelberger, K. L., & Giavalisco, M. 2001, ApJ, 554, 981 
Rauch, M. 1998, ARA&A, 36, 267 

Rauch, M., Haehnelt, M. G., & Steinmetz, M. 1997, ApJ, 481, 601 
Rauch, M., Sargent, W. L. W., & Barlow, T. A. 2001, ApJ, 554, 823 
Schaye, J. 2001, ApJ, 559, 507 

Schaye, J., Rauch, M., Sargent, W. L. W., & Kim, T. 2000a, ApJ, 541, LI 
Schaye, J., Theuns, T, Rauch, M., Efstathiou, G., & Sargent, W. L. W. 2000b, 

MNRAS, 318, 817 . . 

Springel, V. & Hernquist, L., 2002, preprint ( |astro-ph/0206393| ) 
Thacker, R. J „ Scannapieco, E , , & Davis , M . ^OU^, ApJ, submitted (p stro- 



ph/0205054X" 
Theuns, T, Leonard, A.. 



, Efstathiou, G., Pearce, F. R., & Thomas, P. A. 
MNRAS, 301, 478 
Theuns, T, Schaye, J., & Haehnelt, M. G. 2000, MNRAS, 315, 600 
Theuns, T, Mo, H. J., & Schaye, J. 2001, MNRAS, 321, 450 
Theuns, T., Bernardi, M., Frieman, J., Hewett, P., Schaye, J., Sheth, R. K. 
Subbarao, M. 2002, ApJ, 574, L 1 1 1 



1998, 



Theuns et al. 



5 



Vial M, Matarrese S.Mo HJ, Theuns T, Haehnelt MG, 2002, MNRAS in press 
Ustro-ph/0203418) 



Webb, J. K. 1987, PhD thesis, Univ. Cambridge 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 



